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Abstract 

For studying the dynamics of a two-level system coupled to a quantum oscillator we 
have presented an analytical approach, the transformed rotating-wave approximation, 
which takes into account the effect of the counter-rotating terms but still keeps the 
simple mathematical structure of the ordinary rotating-wave approximation. We have 
calculated the energy levels of ground and lower-lying excited states, as well as the 
time-dependent quantum dynamics. It is obvious that the approach is quite simple 
and can be easily extended to more complicated situation. Besides, the results are 
compared with the numerically exact ones to show that for weak and intermediate 
coupling and moderate detuning our analytic calculations are quantitatively in good 
agreement with the exact ones. 
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I Introduction 



The physics of a two-level system coupled to a quantum oscillator (spin-oscillator 
model, SOM) is of wide interest because it provides a simple but ubiquitous model for 
numerous physical processes, such as the superconducting qubit of Josephson junction[l, 
2, 3, 4, 5, 6], the semiconductor quantum dot [7, 8], the coupling between a qubit and 
a nanomechanical oscillator [9, 10], and a toy model for Holstein polaron[ll]. The 
Hamiltonian of SOM reads 

H = ln<j x +u;tfb + ?-(tf + b)a z . (1) 

fl is the level difference and a x and a z are Pauli matrices to describe the two-level 
system, (b) is the creation (annihilation) operator of the quantum oscillator with 
frequency uj and g the coupling constant between the two-level system and the oscil- 
lator. The model seems quite simple, however, an analytical solution has not yet been 
found and various approximate analytical and numerical methods have been used. 

The Hamiltonian (1) is equivalent to the famous Jaynes-Cummings model with 
inclusion of both the rotating- wave terms and the counter-rotating terms [12, 13], 

Hj C = + outfb - 9 -<j x {tf + b), (2) 

if a rotation around the y axis is taken for the Pauli matrices, e l7Tav ^a x e~' t7ray ^ 4 = 
a z and e Mro " !/ / 4 o- 2 e _ " rcr! '/ 4 = — a x . Then, (b) is the creation (annihilation) operator 
of the cavity mode and Q is the atomic transition frequency (g is the vacuum Rabi 
frequency). Becides, the Hamiltonian (1) is also of interest in the research area of 
matter-field interaction because it is the following atom-field interaction Haf in the 
rotating frame[8, 14], 

Haf = \w Q a z + ^ (a+e^ + as™*) + utfb + f (& f + b)a s , (3) 

where a± = ^(a x ±ia y ) . Here uo is the atom transition frequency (or the exciton energy 
in quantum dot) and uj = u is the frequency of laser field (in resonance with cu ), and 
fl the Rabi frequency. The dynamical evolution of the interacting Hamiltonian Haf 
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can be described by 

p(t) = {^AM{t)w z \4>AM{t)) = (mWzim), (4) 

where |t/>am(£)) is the wave function of Ham and \ip(t)) = e Wst(Jz ^\i)AM{t)) is that of H 
(Eq.(l)). The initial state of the system is assumed to be \ip(0)) = \i/)am(Q)) = I T)|0), 
where a z \ f) = ] |) is the eigenstate of a z and |0) is the vacuum state of the quantum 
oscillator. When the coupling g — it is easily to get the typical Rabi oscillation 
P{t) = cos(At). For g ^ 0, the Rabi oscillation may be modulated by the interaction 
with the quantum oscillator and it is an interesting problem related to the quantum 
manipulation of the interacting system[6, 10]. 

Although there is still no analytically exact solution for (1) or (2), various approx- 
imate analytical solutions for (1) or (2) already exist[9, 10, 15, 16, 17, 18]. The most 
popular one may be the rotating- wave approximation (RWA), which relies on the as- 
sumption of near resonance (\Q — uj\ <C and u) and weak coupling (5 < and a;). 
However, the condition of near resonance may not be satisfied if we start from Hamil- 
tonian (3), because Q is the Rabi frequency of the pumping field which may be much 
larger or smaller than the oscillator frequency [8, 14]. Besides, quantum-limited solid- 
state devices are an alternative to the usual atom-cavity implementation of the SOM, 
and in these solid-state devices the coupling strength g and the detuing \Q — uj\ may be 
outside the regime where the RWA is valid (the coupling between a nanomechenical res- 
onator and a charge qubit may be in the intermediate range g/uo m 0.1 ~ 1[6, 10, 18]). 
This is the motivation for approximate analytical solutions beyond the RWA to be 
presented in last years. 

The Hamiltonian (1) or (2) can be numerically solved easily and quickly by ordinary 
PC, then, why do we still need an approximate analytical solution? As far as we can 
see, the purpose may be: l.to see more clearly the physics, 2. to test the accuracy of 
the analytical solution for extending it to more complicated model where a numerically 
exact solution is difficult to be obtained. Hence, we have following criterion for the 
validity of an approximate analytical solution: First, it should be as simple as possible, 
then, it can be easily extended to more complicated situation. Second, the main physics 
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should be catched, at least for the interesting and concerned range of the parameters, 
and it should be as accurate as possible compared with the numerically exact result. 

Recently, one of us proposed an analytic approach[19] to the spin-boson model, 
which describes a two- level system coupled to a dissipative environment [20]. Roughly 
speaking, Hamiltonian (1) is a simplified version of the spin-boson model, that is, the 
single mode spin-boson model. Our approach, which is a perturbation expansion based 
on the unitary transformation, has been successfully applied to several problems related 
to the interaction between the quantum-limited system and its environment [21, 22, 14]. 
In this work, we would show the validity of our analytic approach by studying the nu- 
merically solvable model, which leads to the transformed rotating-wave approximation 
(transformed RWA) for model (1). We will focus on the dynamical evolution P(t), since 
the accuracy of its calculation depends not only on the calculation of ground state but 
also on that of the lower-lying excited states. The dynamical evolution P(t) calculated 
by our transformed RWA will be compared to the numerically exact one, as well as 
to the ordinary RWA, to show that for weak and intermediate coupling and moderate 
detuning our calculations are quantitatively in good agreement with the numerically 
exact results. 



II Theoretical analysis 

The ordinary RWA for the SOM (Eq.(l)) is 



I g p -, 

Hrwa = -^o- x +utfb+ - \w\si)(s 2 \ +b\s 2 ){si\\ , (5) 

where <7 x |si) = — | s i) and u x \s 2 ) = \s 2 ) are the eigenstates of a x . It is easily to see that 
|si)|0) (|0): the vacuum state of the oscillator) is the exact ground state of Hrwa an d 
the mathematical structure of Hrwa is quite simple. 

We apply a unitary transformation[19] to H (Eq.(l)), H' = exp(S)H exp(— S), and 
the purpose of the transformation is to take into account the effect of counter-rotating 
terms, where 

S = b). (6) 
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A parameter £ is introduced in S and its form will be determined later. The transfor- 
mation can be done to the end and the result is 



H' = H' + H' 1 + H' 2 , (7) 
H' = ^Qa x + utfb - ££(2 - 0, (8) 

H[ = \g{\ - i)(b ] + b)a z + l v nia y ^ ~ b), (9) 

Z Z OJ 

H' 2 = ±na x (cosh{fatf-b)}-r?) 

+±m* v (sinh{^(&t _ b )} - V ^(b j ~ b)) , (10) 

Z \ OJ OJ J 



where 



V = eM-£^e]- (11) 

Obviously, H' Q can be solved exactly because for which the spin and the oscillator are 
decoupled. If the displacement parameter £ is determined as 

( = (12) 
oj + i]U 

then we have 

H[ = \b\a z + ia v ) + b(a z - ia y )} = r]Q^ [^siX^I + &|s 2 )(si|l . (13) 



2 oj 



OJ 



Note that H[ is of the same form as the RWA in Eq.(5), except the different coefficient. 

The transformed Hamiltonian H' is equivalent to the original H and there is no 
approximation till this point. In the following, the transformed Hamiltonian is approx- 
imated as H' ks H TRWA = H' + H[, since the terms in H 2 are related to the double- 
and multiple-boson transition and their contributions to the physical quantities are 
0(g 4 ). That means, through the unitary transformation approach we get the trans- 
formed RWA Hamiltonian Htrwa, which is of the same mathematical structure as the 
ordinary RWA H RWA in (5). |si)|0) is the exact ground state of H TRW a with ground 



state energy 



tf„ = -^n-|^(2-o. (14) 



The eigenenergy for excited states can be easily obtained from H[ (Eq.(13)), since 
it contains the rotating- wave terms only. For n = 0,1,2,..., the eigenenegies for all 
excited states are 

E 2n+ i = (n + - i^-^ + ^n + l) - £^(2 - 0, (15) 
E 2n+ 2 = (n + \)u + y^-rfflt + gvin+l) - £^(2 - £), (16) 

where g' = 2gi]Q/(uj + rjQ). 

The calculations of the ground state and lower-lying excited states with moderate 
detuning are shown in Figs.l and 2. Fig.l is for the case of Q = 1 and uj = 0.5 but 
Fig. 2 for = 0.5 and u — 1. For comparison, the results of exact diagonalization and 
those of generalized RWA[10] are also shown. For ground state energy our result of 
transformed RWA is very close to the exact result and is better than that of generalized 
RWA. For excited states we could not say definitely which one is better, but we can 
say that for intermediate coupling g/u & 0.1 ~ 1 our results of transformed RWA are 
quite close to the exact results. 

Ill Dynamical evolution P(t) 

We have shown that for the weak and intermediate coupling our calculation of the 
ground state and the lower-lying excited state is quantitatively in good agreement 
with the numerically exact results, which is a check of the validity of our approach. 
Furthermore, in this section our approach will be checked by calculation of the excited 
state properties, that is, the dynamical evolution P(t) of Eq.(4). 

The numerical exact calculation for P(t) with Hamiltonian H can be obtained by the 
numerical diagonalization. The following is the calculation of P(t) in our transformed 
RWA. From Eq.(4) we have 

P(t) = (iP(0)\e iHt a z e- iHt \iP(0)) = (^(0)|e" s r e iH V : a z e~ s ' e -^V|V>(0)) 
w (ip'{0)\exp(iH TRWA t)(T z exp(-iH TRWA t)\ilj'{0)) 

= (^(OI^W-^KW), (17) 
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where \ip'(0)) = e s \^(0)). Here the unitary transformations, e s He s = H' and 



e a,e 



-s 



a z , have been used. IV'j(O) = etH ° l e w{—iH T RWAt)W{fy) is the wave 
function in interaction picture, which is the solution of following Schroedinger equa- 
tion, 



9' 



(18) 



H[(t) = e iH ^H[e- iH ^ = | (tf\ Sl )(s 2 \e l( ^ n)t + b\s 2 ) (s^e'^-^) . (19) 

The initial state is |^'(0)) = e s | |)|0) = + |s 2 ))e Q{bt " 6) |0) with a = g£/2uj. 

Here the exponential operator is expanded to the order a: e a ^' h ^ « 1 + a(b jf — b) and, 
hence, the initial condition is: 

Cio(O) = C 20 (0) = ± C n (0) = C 21 (0) = C 12 (0) = 0, (20) 

with the wave function [23] 

|^(*)) = C 10 (t)|s 1 )|0) + C7 2 o(*)|s2>|0> 

+Cn(t)|si>|l>+C7 2 i(*)|s2)|l> + Ci 2 (*)|si)|2). (21) 



The Schroedinger equation is [23] 



»4<3»W = ^e>W-"»C n (t), 4c„(i! V -""- , -"<:\, 1 (f). 
at 2 at 

(I .. q 1 ,,„o . d 



9_ p -i( v n-u>)t ( 
2 



*^ 2l(t) = f^-^Ci^), i^C 12 (f) = | e -*^-")«C7 21 (0. 



Then, the solution for Eq.(23) is 



<3> t i]Q -u $ t 

cos z — sm 

2 $ 2 



— a — sm > e 

$ 2 J 



i(r)Q-u>)t/2 



§ t r)tt-u § t 

cos h « — sm 

2 $ 2 



-^-sin^U-^)*/ 2 , 
$ 2 J 



and the solution for Eq.(24) is 



C 21 (t) 



a 

72 



cos ? — sm 

2 $i 2 



c 12 (t) = 



a i 



V2g' . $it 



sin e 

y/2 $i 2 



(22) 
(23) 
(24) 

(25) 
(26) 

(27) 
(28) 



where $g = (rjQ - u) 2 + g' 2 and = (77ft - cj) 2 + 2g' 2 . 
The dynamical evolution is 

P(t) = 2Re (C 2 * (i)C 10 (Oe^ + C^C^e^) . (29) 

For the ordinary RWA (Eq.(5)), the interaction is 

Hl (t) = | (6t| Sl )( S2 | e ^^ + 6|s 2 )( Sl |e-^-^) . (30) 

Then, it is easily to get the dynamical evolution, 

P R wA{t) = cos cos — sin sin , (31) 

where ^\ RWA = (ft - to) 2 + g 2 . 

The calculations of P{t) with moderate detuning and intermediate coupling are 
shown in Figs. 3 and 4. Fig.3 is for the larger ft case (ft = 1, u — 0.5 and g = 0.4) but 
Fig. 4 for larger lo case (ft = 0.5, u — 1, and g = 0.5). For comparison, the results of 
exact diagonalization and those of ordinary RWA are also shown, one can see that our 
results of transformed RWA are quite close to the numerically exact ones, much better 
than those of ordinary RWA (especially for the case of Fig.4). 

For some special case the ordinary RWA is absolutely invalid. Fig. 5 shows the P(t) 
dynamics for large detuning lo — 1, ft = 0.25 and strong conpling g = 1. One can see 
that the result of ordinary RWA is qualitatively incorrect, but that of our transformed 
RWA is still in qualitative agreement with the numerically exact one. 

IV Concluding remarks 

We have presented an analytical approach, the transformed RWA, for the SOM, which 
takes into account the effect of counter-rotating terms but still keeps the simple math- 
ematical structure of the ordinary RWA. We have investigated the energy levels of 
ground and lower-lying excited states, as well as the time-dependent quantum dynam- 
ics. It is obvious that the approach is quite simple and can be easily extended to 
more complicated situation. Besides, the results are compared with the numerically 
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exact ones to show that for weak and intermediate coupling and moderate detuning 
our analytic calculations are quantitatively in good agreement with the exact ones. 

Generally speaking, the generator (6) is for a displacement transformation and the 
parameter £ may be used to indicate the amount of the displacement of the oscillator. 
If £ = 0, that is, without the transformation, we have the ordinary RWA (Eq.(5)) which 
is good for weak coupling and near resonance. If £ = 1, that is, the oscillator can follow 
adiabatically the tunnelling of the two-level system, we can have the generalized RWA 
of Ref. [10] which is good for higher frequency uj Q and strong coupling g/uj 1. 
Our £ is determined in Eq.(12), < £ < 1, which is in between the ordinary RWA and 
the generalized RWA. Physically, our £ is to take into account the nonadiabatic effect 
when the retardation of the interaction between the two-level system and the quantum 
oscillator with moderate detuning and intermediate coupling is important. In addition, 
the bare coupling g in original Hamiltonian H is replaced by the renormalized coupling 
g' = 2gr]Q/(u + i]Q) in Htrwa because the effect of the counter-rotating terms has 
been included. 
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Figure Captions 



Fig.l The energy levels of the ground state (E g ) and the lower-lying excited states 
(Ei, E 2 and E 3 , from bottom to top) as functions of the ratio g/uj. Q = 1 and uj = 0.5. 
Our results (dashed lines) are compared with the numerically exact ones (solid lines) 
and those of the GRWA Ref.[9] (dashed-dotted lines). 

Fig. 2 The energy levels of the ground state (E g ) and the lower-lying excited states 
(Ei, E 2 and E 3 , from bottom to top) as functions of the ratio g/uo. VL = 0.5 and uj — 1. 
Our results (dashed lines) are compared with the numerically exact ones (solid lines) 
and those of the GRWA Ref.[9] (dashed-dotted lines). 

Fig. 3 Time evolution P(t) for Q = 1, u — 0.5 and g = 0.4. The dashed line is 
our transformed RWA calculation. The solid line and the dashed-dotted line are the 
numerically exact result and the ordinary RWA one, respectively. 

Fig. 4 Time evolution P(t) for VL = 0.5, uj — 1 and g = 0.5. The dashed line is 
our transformed RWA calculation. The solid line and the dashed-dotted line are the 
numerically exact result and the ordinary RWA one, respectively. 

Fig. 5 Time evolution P(t) for f2 = 0.25, uj = 1 and g — 1. The dashed line is 
our transformed RWA calculation. The solid line and the dashed-dotted line are the 
numerically exact result and the ordinary RWA one, respectively. 
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